1 



QCD at finite baryon density with t-asymmetric fermions * 

E. Mendel, L. Polley ^ 

^FB Physik, Universitat Oldenburg, 
D-26111 Oldenburg, Germany 

Susskind's continuous-time fermions, with two flavours, can be latticized using a one-sided time derivative. We 
are presently investigating the interacting case, where we hope to find the onset at finite fi at the right place due 
to the reduced number of flavours. As for these fermions there is only a discrete chiral symmetry left over, the 
lightness of pious in the broken phase has to be investigated. 



1. INTRODUCTION 

Based on the idea that the early onset of 
Baryon density as function of chemical poten- 
tial, /J,, could be due to the excessive number of 
valence-quark flavours present in previous simula- 
tions, we have started to investigate t- asymmetric 
fermions, where we can reduce to two flavours. 

In the past, several [1-3] lattice QCD simu- 
lations at finite jj, had shown an onset in ther- 
modynamic quantities at a much lower than 
at the expected mf^/S. One of us had proposed 
[4,2] that this early onset could be due to the fact 
that we had been using the usual Kogut-Susskind 
fermions, thus having a 4-fiavour valence-quark 
degeneracy, in contrast to the 2 light fiavours 
present in nature. The lattice fermions can 
change their fiavour content along the propaga- 
tion in the interacting case, producing in this way 
many more light nucleon states than in nature, 
thus weakening the Fermi repulsion, while on the 
other hand increasing the scalar Yukawa attrac- 
tion due to the extra light scalar channels. Both 
effects would increase substantially the nuclear 
binding, producing the observed early onset on 
the lattice. 

In order to reduce the quark fiavour degener- 
acy we first suggested the use of Wilson fermions, 
where one has just 1 fiavour, but they are harder 
to implement due to the Dirac components be- 
sides their explicit breaking of chiral symmetry. 
That project is being done with an innovative vol- 
ume method to obtain the fermionic expectation 
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values and its first results are reported in this 
proceedings [5]. 

Here we would like to present another fermion 
formalism, which reduces to 2 the fiavour number 
[6], obtained by using a t-asymmetric fermionic 
action (one-sided time derivative). This action 
is unitarily equivalent to Susskind's Hamiltonian 
formalism [7,8], but in a simpler notation. The 
Heidelberg collaboration [9] were the first to in- 
vestigate a finite /j, for these fermions, in a version 
with 4-component spinors. 

We have shown that at least in the free case 
[6] these fermions behave as 2 genuine fiavours. 
We have constructed their Fock space and proven 
that in the continuum limit they will contain par- 
ticles and antiparticles with the usual spectrum 
for the Hamiltonian and the Baryon number op- 
erators. Furthermore, the transfer matrix can be 
made positive, in contrast to the usual staggered 
fermions case. 

In this work, we will introduce the formalism 
in the SU (3) interactive case and present the first 
simulation results. 

2. t-ASYMMETRIC FERMIONS 

Usually we discretize the fermionic action, Sp, 
with a symmetrical derivative 

^"^ = ^ + i>) - ipix - £>)) 

in order to get good hermiticity properties for Sp- 
This produces the fiavour doubling in each coor- 
dinate direction. But in order to have a sensi- 
cal theory we only need an hermitian hamiltonian 
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(eventually, on the lattice, a positive transfer ma- 
trix). This means that we don't have to take a 
symmetrical time derivative, thus reducing by 2 
the flavour number. 

We choose then an action Sp = '4'xKxy'4'y j 
where: 

K^y = m6^y + 74 (e^ U^^4 6^^y_4- 6^,y^ (1) 

k 

To reduce further the flavour degeneracy, we 
can absorb a 74 into the ip functional integra- 
tion and then, diagonalizing spinors in analogy 
to Kawamoto-Smit: 

-ip(x) = al^af/al^x{x) with at = ij4lk (2) 
we get, 

Ky = mj4T4x)6^y + (e''?7^,4 6^^y_^ - 6^^y^ 

-2t W (^-^ ^.,y-i^ - ^Uk ^...+0 (3) 

k 

with the usual: T^{x) = -K^^-i . Taking 

just the flrst component of this diagonal kernel 
K' (like the thinning out for K.S. fermions) we 
have managed to reduce to just 2 quark flavours. 
Note that our fermions are unitarily equivalent 
to Susskind's hamiltonian formalism ones [7], by: 

Xour = i"+^+"(-ir"Xtheir. 

Susskind showed [7] that there is a discrete 
version of chiral invariance left over, preventing 
renormalization mass counterterms in the inter- 
acting case. This is a nice feature as we know then 
that the chiral limit in the continuum limit should 
be at niq = 0, in contrast to Wilson fermions. Un- 
fortunately, although the chiral condensate seems 
to be spontaneously broken from our flrst simula- 
tions (also at strong coupling [8]), as the chiral 
symmetry is only discrete, there is no guaran- 
tee that the tt will be a Goldstone boson. We 
see that these fermions have some "intermedi- 
ate behaviour" in its properties, in between usual 
Kogut-Susskind and Wilson ones. 

Even though there is no usual Goldstone mech- 
anism, we could still get light pions, in a sim- 
ilar fashion to Aoki's proposal [10] for Wilson 



fermions: Close to a second-order phase transi- 
tion, to a phase where the order parameter (for 
parity), (tt), gets an expectation value, the cor- 
relation length for (7r(0)7r(t)) diverges, giving a 
massles pion. In the Wilson case one can shift k 
in order to get close to the chiral limit, but here 
we have to get the critical point at niq = 0, which 
will presumably only be reached in the continuum 
limit. 

A potential problem with these t-asymm. 
fermions, in contrast to the usual formulations, 
lies in the fact that the action Sp is not mani- 
festly 0(4) invariant, endangering relativistic co- 
variance. One could even naively ask from a 
hopping parameter expansion, if there is back- 
ward propagation at all! In fact, for free fermions 
we have checked that the antiperiodic boundary 
condition plays a crucial role in making positive- 
energy modes propagate forwards and negative 
backwards, in a symmetrical way in the contin- 
uum limit (for modes far from the cutoff). It 
will be important to check that this also happens 
when calculating a meson propagator as it con- 
sists of a particle and an antiparticle propagator. 

Furthermore, even if the discretization in the 
action looks naively symmetric for spatial and 
time directions (when = a^), the fermion 
propagator has to be renormalized (even in the 
quenched case by gluon loops) which will almost 
certainly destroy the 0(4) symmetry. Assum- 
ing that the action is renormalizable (which is 
as yet not clear), we would get renormalization 
constants and Zr in front of the spatial and 
temporal derivative terms in Sp ■ One can absorb 

in the y flelds but then one gets a factor 
in front of the spatial term. Thus, in order to re- 
gain 0(4) invariance we could tune the parameter 
^ = as if we were on naively asymmetric lat- 
tices, until we regain the symmetry. This could 
be controlled non-perturbatively by calculating 
some physical propagator in several directions or 
by evaluating for some states the dispersion rela- 
tion: 




'^{pja„) + {marf (4) 



so as to approach the relativistic one for some ^. 
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3. STATUS OF COMPUTATIONS 



Baryon Density vs Chemical Potential 



We have already done first simulations of the 
baryon density, (J4) = {ef Xx U:c,aXx+'a) the 
chiral condensate, (xx); function of ji. The 
operators have been calculated with the solved 
pseudofermion method [2]. The curves for (J) 
are shown in fig. 1 and show a clear onset for both 
niq for jj, 0.2. This is in clear contrast with our 
results for staggered fermions where there was a 
clear shift among the onsets for both masses. If 
the onset would be controlled by m^/2 (as argued 
by some authors [11]), and if the pton is a Gold- 
stone boson for these fermions, then the onset 
should be at twice /j, for the higher mass. This 
shows that it is crucial to calculate the hadron 
masses in order to decide if the onset is at mjv/3 
as wanted, or at some lower value due to some 
non-Goldstone meson. The vacuum value [6] of 
(J) = 1.5 corresponds to saturation of the Fermi 
sea for these fermions ( 2j^ x 2^ x 3c per cube of 
size 2^). 

The hadron masses can be calculated from the 
decay of the corresponding zero-momentum cor- 
relation of operators [8]. For the tt: 

T^oir) = xt Ur,z Xr+z + xj+i t^r+i,^ Xr (5) 

and for the neutron with spin down, e.g., 

= Hi - i-y + i-T^' + X (6) 

[1 - i + (-)^(l + i)] {-T^'-eABcX^x'^x; 

We are presently starting to calculate these 
masses and studying the 0(4) symmetry on im- 
proved algorithms in order to get the necessary 
precisions. 

4. CONCLUSIONS 

We have seen that for both 's the onset is at 
very similar //'s indicating that it could be con- 
trolled by mjv/3 as wanted. It could also be that 
the TT is not Goldstone-like for these fermions, 
thus being able to give so similar onsets. We are 
currently investigating the masses and renormal- 
ization effects in order to decide this question. 
We would like to thank the RRZN in Hannover 
for supercomputer time. 
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Figure 1. (J) vs /j, for two rriq. The nearby onset 
is inconsistent with a Goldstone-like m^/2. 
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